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We study the single transverse-spin asymmetry for dijet production in hadronic cohisions in 



. both the cohinear QCD factorization approach and the Brodsky-Hwang-Schmidt model. We show 

^: 

JL , that a nonvanishing asymmetry is generated by both initial-state and final-state interactions, and 
M-i, 

(— I ' that the final-state interactions dominate. We find that in the leading kinematic region where 



the transverse momentum imbalance of the two jets, q± = Pi± + P2±, is much less than the 
momentum of either jet, the contribution from the lowest non-trivial perturbative order to both 
the spin-averaged and the spin-dependent dijet cross sections can be factorized into a hard part 
that is a function only of the averaged jet momentum P± = {Pi± — -P2±)/2, and perturbatively 
. generated transverse momentum dependent (TMD) parton distributions. We show that the spin 

asymmetry at this non-trivial perturbative order can be described by the TMD parton distributions 
defined in either semi-inclusive DIS or the Drell-Yan process. We derive the same hard parts from 



■ both the collinear factorization approach and in the context of the Brodsky-Hwang-Schmidt model, 

verifying that they are not sensitive to details of the factorized long distance physics. 
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I. INTRODUCTION 



Single-transverse spin asymmetries (SSAs) in high-energy colhsions with one transversely 
polarized hadron are important phenomena that have been observed for more than three 
decades in various physical processes Q, 13, [s, Q, [s, [g], [t]. In these processes, the observed 
final-state hadrons show an asymmetric distribution in a plane perpendicular to the beam 
direction when the transversely polarized hadron scatters off an unpolarized hadron (or a 
virtual photon). The SSA is defined as An = {<7{S±) — cr{—S±))/{(T{S±) + cr{—S±)), the 
ratio of the difference and the sum of (differential) cross sections when the hadron's spin 
vector, S±, is flipped. Recent experimental measurements of SSAs in polarized semi- inclusive 
lepton-nucleon deep inelastic scattering (SIDIS) 3, 4| and in hadronic collisions 0, [g], have 
renewed the theoretical interest in SSAs and in understanding their roles in hadron structure 
and Quantum Chromodynamics (QCD). Although it was realized a long time ago [8;] that 
perturbative QCD can be used to study the effects of transverse spin, the size of the observed 
asymmetries came as a surprise and has posed a challenge for researchers in this fleld jsj. 

Within a model calculation [l^ , Brodsky, Hwang, and Schmidt showed that the flnal state 
interaction in deep inelastic scattering (DIS) can generate a phase required for a nonzero SSA 
in SIDIS. It was later realized that this flnal state interaction can be factorized into the gauge 
link of the gauge invariant transverse momentum dependent (TMD) quark distributions 
The nonvanishing SSA obtained in Ref. [l^ is a consequence of the existence of a naively 
time-reversal-odd TMD quark distribution, the so-called Sivers function . When applying 
the same calculation to the Drell-Yan production of lepton pairs at hadron colliders, the flnal 
state interaction in SIDIS becomes an initial state interaction in Drell-Yan, and the phase 
changes sign, which leads to a prediction of a sign change in the SSAs between these two 



processes m^- This uoutriva. "uuiversahty" property associated witlr the TMD partou 
drstnbutions . the co_ce of gauge Weractions rn QCD (li y, M, and of the QCD 

factorization theorems for these two processes |15l . 1161 . 1171 . 118I | . Experimental tests of this 



prediction will be crucial for our understanding of the origin of SSAs in QCD 19|. 

In Ref. {2^, it was proposed to study the Sivers functions by means of a SSA in azimuthal- 
angular correlations of two jets produced nearly back-to-back at hadron colliders. Mea- 



surements of this SSA for dijet production have begun at RHIC 



2l| . complementing the 



measurements in SIDIS. Unlike the SIDIS or Drell-Yan process, dijet production at hadron 
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colliders involves both initial and final state interactions that may produce the phase needed 
for a nonvanishing SSA. Consequently, the sign and the size of the asymmetry will depend 
on the relative strength of these interactions. Following the previous works on SIDIS and 
Drell-Yan, the authors of Ref. {22! developed a systematic approach to describe the role of 
initial- and final-state interactions in generating SSAs in hadronic collisions, and they found 
that summing all initial/final state interactions into the gauge link of the TMD parton 
distributions leads to a very complicated functional form of the gauge link. In particular, 
the TMD parton distributions studied in dijet correlations in hadronic scattering will have 
no connection to those in the SIDIS and Drell-Yan processes, because their definitions are 



completely different [22|. One then has to question the universality of the TMD parton 
distributions, and the predictive power of perturbative QCD calculations, which relies on 
comparing physical observables with the same factorized long-distance physics while having 
different perturbatively calculable short- distance dynamics. 



In Ref. 231], we briefly reported a new result for the SSA in dijet production in hadronic 



collisions in the twist-3 Efremov-Teryaev-Qiu-Sterman (ETQS) approach 2^, |25|]. We con 



sidered the spin-dependent cross section, Aa{S±) = {cr{S±) — (t{S±))/2, for the process 

A{Pa, S^) + B{Pb) ^ Ji(Pi) + J2{P2) + X , (1) 

with the jet momenta Pi = P + q/2 and P2 = — P + q/2. When both Pi and g_L, the 
transverse components of the momenta P and g, respectively, are much larger than Aqcd? a 
nonvanishing A(t(S'i) can be generated by the ETQS mechanism in the coUinear factoriza- 
tion approach. In this framework, the SSAs are attributed to the spin-dependent twist-three 
quark-gluon correlation functions, which correspond to a quantum interference between dif- 
ferent partonic scattering amplitudes. Since the incoming partons are approximated to be 
collinear to the corresponding incoming hadrons in this approach, the momentum imbal- 



ance of the two jets is generated by producing a three-parton final-state. In Ref. [23|], we 
calculated the contribution from initial-state gluon radiation to Acr(5'i) in the kinematic 
region where Pi ^ gj. ^ Aqcd- We presented the final result for the leading contributions 
to Acr(S'i) in the expansion of the partonic scattering in q±/P± involving a hard qq' qq' 
subprocess. In this paper, we will provide the detailed derivations of this result, and we will 
also present the full contributions from all other partonic subprocesses at the same order. 
We find that although both initial-state and final-state interaction lead to a nonvanishing 
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SSA, the final-state interactions give the dominant contributions to Acr(5'_L). We therefore 
expect that the SSA in dijet production will have the same sign as the Sivers asymmetry in 
SIDIS. 

We find that at leading order in the q^jP^ expansion, the perturbatively calculated 
partonic parts can be further factorized into a single-scale (-Pl) hard part and perturbatively 
generated TMD parton distributions with transverse momenta = (9(g_L). We also find 
that our perturbatively calculated result is equal to the leading-order term in the Aqcd/q'j 
expansion of the following generalized TMD factorization formula 



23|: 



ab 



h± ■ q± 

Mp 



SIDIS / 



(2) 



where Ylia b ^^^^ '^^^^ parton flavors, H^l^^cd Sabred are partonic hard and soft factors, 
respectively, and the [ ]c represents a trace in color space between the hard and soft factors 
due to the color flow into the jets 26|, |27|]. The hard factor in Eq. ([2]) only depends on the 
single hard scale Pi in terms of partonic Mandelstam variables of the reaction ab cd: 



with Xa 



oy2^ 



u 



Xb 



{Pa+Pbf 
{Pa - Pcf 
{Pb - Pcf 



-Pl [ey^-y^ + 1) 



,2/1-1/2 



1) 



(3) 



(e + e and yi and y2 the rapidities of the two jets. 



In Eq. ([2]), q^^^^ and ff^^^^ denote the transverse-spin dependent TMD quark distributions 
(known as the Sivers function) and the unpolarized TMD quark distribution, respectively; 
these TMD parton distributions were chosen to follow their definitions in the semi-inclusive 
DIS process. For example, for a polarized proton with momentum P = (P+,0~,0_l) with 
P± = l/v/2 (P° ± P3) and transverse spin vector S±, the TMD distribution f or q uark flavor 
a can be defined through the decomposition of the following matrix element 281], 

p+rfr d% 



Ma 



IT 



{2n) 



{PS\^|J,{0^ioo■,OUoo■,0)MO)\PS) 



(4) 



where Mp is a hadronic mass scale introduced to keep the TMD distributions fa and qxa at 
the same mass dimension, and the gauge hnk C is defined in a covariant gauge as £t,(oo; S,) = 
exp (^—ig dXv ■ A{\v + ^)) with the path hnk extended to +00. t> is a vector conjugate 
to the momentum vector P. Since we wiU work in a covariant gauge throughout this paper, 
the vector v could be chosen to be a light-cone vector with v"^ = and v ■ P = 1. An 
off-light-cone vector (f ^ 7^ 0) will have to be used when high order corrections are taken 



into account 17|. If we work in a singular gauge, like the light-cone gauge, an additional 
gauge link at the spatial infinity = +00) will have to be included in order to ensure the 
gauge invariance of the above definitions |l3|. We have chosen TMD parton distributions 
defined in SIDIS because of the dominance of final-state interactions. Choosing TMD parton 
distributions defined according to the Drell-Yan process would change the sign of the partonic 
hard factors, but not affect the overall sign of the physical cross section. We note that the 
various factors in Eq. (jlj), apart from being functions of transverse momentum, also depend 
on the renormalization and factorization scales, and especially on the gluon rapidity cut-off 
when higher-order corrections are taken into account 15, Q, The latter dependence 



which leads to a resummation of 



is governed by the Collins-Soper evolution equation 

large logarithms of the form a"ln^"~^ [P±/<1±] in the perturbative series |l5l. Il6l|. 

In order to evaluate the SSA for dijet production, we also calculate in this paper the 
contributions from initial-state gluon radiation to the spin-averaged dijet cross section, a = 
{cr{S±) + a{S±))/2, in the same kinematic region where P± ^ q±. Using the coUinear 
factorization approach, we find that like the spin-dependent case, the leading contribution 
in the q±/P± expansion of the perturbatively calculated partonic scatterings can be further 
factorized into a hard part at a single-scale Pi and perturbatively generated unpolarized 
TMD parton distributions with transverse momenta k± = 0{q±). Our leading contribution 



in the q±lPx expansion is equal to the leading-order term of the following TMD factorization 
f 1 
291, 



formula 



dyidy2dPld^q± ^ J 

[Sat^M±)H:;^^cAPl)], S^'Hka± + h± + X±- q±) , (5) 



ah 
X 



where the superscript "um" indicates the scattering of an unpolarized beam off an unpolar- 
ized target/beam. For consistency, we again express the factorization formula in terms of 
TMD parton distributions defined in the SIDIS process. Actually, the unpolarized TMD 



parton distributions defined in SIDIS and the Drell-Yan process are identical, because the 
unpolarized parton distributions are invariant under the naive-time-reversal transformation. 

A key feature of any QCD factorization is that the perturbatively calculated short- 
distance hard factors should not be sensitive to details of the factorized long distance physics. 
In this paper, we will derive all short- distance hard factors in Eqs. 02]),® also by using the 
Brodsky- Hwang- Schmidt model for SSAs [l^. We find that the hard-scattering factors 
derived in this model are in fact the same as those derived in the coUinear factorization ap- 
proach, despite of the clear difference in the treatment of the nucleon. Both the spin-averaged 
and the spin-dependent cross sections calculated in this model are therefore consistent with 
the generalized factorization formulas in Eqs. 

To fully investigate the above factorization formalism, we would have to consider the 
contributions from gluon interactions and radiations in all possible regions of the phase 
space, and to all orders in perturbation theory, which is not what we are trying to do in 
this paper. Instead, as a step to test the above factorization formalism, we will study the 
contribution at the first non-trivial order with a single gluon radiated nearly parallel to one 
of the incident nucleons. From this study, we will be able to show how to factorize this gluon 
into the TMD parton distributions of the incoming nucleons, and to verify the definition 
of the TMD parton distributions used in the generalized factorization formulas. Our first 
order calculation at P± ^ q± ^ ^qcd clearly shows the factorization of a hard part at 
0{P±) and that the radiation of a collinear gluon can be absorbed into the relevant TMD 
parton distributions. However, our work does not address the factorization between the 
various TMD parton distributions at the scale 0{q±). Our calculation should be regarded 
merely as a "necessary condition" for such a factorization to hold. Since the predictive power 
of perturbative-QCD calculations involving measured hadrons relies on factorization, a full 
proof or disproof of this factorization remains an important challenge [s^. 

The rest of this paper is organized as follows. In Sec. II, we calculate the partonic 
hard factors in a simple model-inspired approach. Using the quark-diquark model for the 
proton and the Brodsky- Hwang- Schmidt model for SSAs, we extract the hard factors by 
the factorization of the TMD parton distributions. We find that the hard factors for the 
unpolarized scattering in Eq. ([2]) are identical to the partonic differential cross sections 
da/di, and the ones for the transverse-spin dependent cross section are different due to the 
initial- and final-state color interactions. In Sec. HI, we calculate the SSA in the collinear 
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a, /i 



FIG. 1: The vertices in the quark- diquark model that we will use to calculate the hard factors. Here 
"a" represents the color-index for the gauge boson, and i,j the color-indices for the quark and/or 
the di- quark. 



factorization approach, and show that the colhnear gluons parallel to the polarized proton 
can be factorized into the Sivers function, while those parallel to the unpolarized proton can 
be factorized into the unpolarized TMD parton distribution. We also demonstrate that the 
hard factors calculated in both approaches are the same and independent of the details of 
the factorized long-distance physics. In Sec. IV, we extend our formulas to the g^-weighted 
SSA, and compare with previous results in the literature. Finally, we summarize our paper 
in Section V. 



II. CALCULATION OF THE HARD FACTORS 

In this section, we will calculate the hard factors at 0{P±) in the generalized factorization 
formulas in Eqs. ([2]), ([5]), based on a model in which an energetic parton scatters off a nucleon 
that is made of a quark and a scalar diquark [10] . We will follow the approach of Brodsky 
et al. and incorporate the correct color interactions between the gauge boson and the 
scalar diquark in the proton wave function. We extract the hard factor by factorizing 
the model-sensitive long-distance physics from the differential cross section into the parton 
distributions. Because the model-dependence was factorized into the parton distributions, 
the extracted hard factors resulting from our calculations do not depend on how the nucleon 
couples to the partons, and therefore, should be model-independent. 
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A. Hard factors for the unpolarized cross section 



In this subsection, we derive the hard factors H^^_^^^ for the unpolarized differential 
cross section in Eq. We will first give a simple example of how the hard factors may 
be extracted in the context of the quark-diquark model of [l^. We will also give a more 
general approach based on power counting techniques to calculate the hard factors directly 
from partonic diagrams. 

To specify the quark-diquark model of the proton, we give in Fig. [1] graphical rules for 
the vertex coupling the quark-diquark to the proton, and the vertex coupling the diquark 
to a gauge boson. In order to correctly take into account the color degree of freedom of the 
strong interactions, we included proper color factors for these vertices. Using this quark- 
diquark model to extract the hard factors not only verifies the fact that the hard factors are 
independent of the dynamics at the scale of the proton, but also provides a relatively simpler 
demonstration of how to factorize the cross sections into the TMD quark distributions and 
corresponding hard factors at the leading order in the q±/P± expansion. 

In this particular model, we calculate the cross section for an energetic parton of mo- 
mentum Pb = (0+, Pg , 0±) scattering off a parton of momentum P4 = (P|, 0, 0^) from the 
(modeled) nucleon, producing two jets of momenta Pi and P2 plus an unobserved particle 
with momentum k' in the final state, as shown in Fig. [21 where the Feynman diagram rep- 
resents the contribution from a quark-quark scattering channel: qq' qq' with the quark 
q from the nucleon. We extract the hard factors by comparing the calculated cross section 
with the factorized formulas in either Eq. ([2]) or ([S]). The diquark in the final state has a 
transverse momentum k'j_ = —{Pi± +P2±) to balance the small momentum imbalance of the 
two jets, q± = Pi_L + P2±. 

According to the factorization formulas in Eqs. ([2]), ([5]), the transverse momentum imbal- 
ance q± of the two jets is a consequence of adding the transverse momenta of the two active 
partons and from the soft factor. Since the incoming quark in this model calculation has 
only a longitudinal momentum component, its transverse as well as longitudinal momentum 
distributions should be given by delta functions, 

fb{xi>,PB±) = S{x, - 1)6^^\Pb±) . (6) 

Similarly, at this order, the soft factor is also a delta function of the transverse momentum, 
i.e., S{X±) = 5^'^\X±). Therefore, the only contribution to the transverse momentum imbal- 
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FIG. 2: T/ie Feynman diagram contributing to the dijet- correlation in the qq' qq' channel, for 
the unpolarized cross section. 

ance of the two jets comes from the transverse momentum of the quark emerging from the 
nucleon, and the generahzed factorization formulas can be reduced to 



dyidy2dPl(Pqx 
for the unpolarized cross section, and 



(7) 



(8) 



dyidy2dPld'^q± Mp 
for the single transverse-spin dependent cross section. 

In order to extract the hard factor from the differential cross sections in Eqs. ([7]), ([8]), 
we will need to calculate the TMD quark (unpolarized and Sivers) distributions within the 



same model. They are given by the diagrams shown in Fig. [3l and are available 
can summarize their results as follows: 

1 — X 



10 



13|. We 



f{x, /cj 



9% 

167r3 M{x, kl 



167r3 



[kl + {xMp + mqf] 



2- ^l-x){my + xMp)Mp^^A^{x,kl) 



(9) 
(10) 



A2(x,A;2) kl A2(x,0) 

where Cp = {N^-l)/2Nc and A'^{x,kl) = kl+ xXl + {1 - x)ml - x{l - x)M^, and Mp, rUg 
and \g are masses for the proton, the quark and the gauge boson, respectively. In the above 
formulas, we have included the color-factors arising from the vertices described in Fig. [T] 

We now present details of the calculation of the contribution of the quark-quark scattering 
channel qq' — > qq' to the unpolarized scattering cross section. We also provide full results 
for all other channels, which can be derived similarly. The Feynman diagram for qq' qq' 



(a) 



(6) 



FIG. 3: The leading order unpolarized quark distribution (a) and Sivers function (h), as given in 
the quark- diquark model \l^ J- 

has been shown in Fig. [2l In the limit q± <^ Pi, the two final-state jets have approximately 
the same transverse momentum, i.e., Pi± ~ P2± ~ P±- We define the following kinematic 
variables: 



Xa2 



3 2/1 



,2/2 



Xbl 



Xh2 



Pi 



-2/1 



-2/2 



(11) 

(12) 



where s = {Pa + PbY- From these, we immediately find the incident partons' momentum 
fractions Xa = Xai + Xa2 and = x^i + Xf,2, and the scalar diquark's momentum can be 
written as 

k' = {l- Xa)PX + (1 - Xb)PB + fci • (13) 



In the evaluation of the cross sections, we use the power counting analysis [3l|, l32], keep only 
the leading power contributions and neglect all higher order corrections in k/Pl, where k 
represents any lower mass scale like Mp, ruq and A^. In this limit, the differential cross 
section for Fig. [2] can be written as 

.2 / . \ 2 



1 



\Mmk'y-xt) 



(14) 



dyidy2dPld'^q± 2s ^167?^^ 

where Ai is the scattering amplitude for the diagram. Expanding the above delta function 
for the phase space integral of k', we find that in the leading power contributions 33| . 

1) , 6{Xa-l) 



1 / 5{xb 

4(1- 



Xn 



Xb 



S{xa — 'i-)S{xb — 1) In 



(15) 



where the "plus" distribution follows the usual definition 3J]. This delta function will help 
to simplify our calculations for contributions from different kinematic regions of k'. For 
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example, if k' is parallel to Pa, which means (1 — Xa) 7^ 0, we only have a contribution from 
the first term of the above expansion, which also implies x;, = 1 in this limit. Furthermore, 
the quark propagator in Fig. [2] reads 

(Pa - k'Y - ml = ~ k^{xa,ql) ' ^^^^ 

where A^(xa, q\) is defined above after Eq. ( |TOl) . and is order of . Combining this propa- 
gator with the delta function expansion, we find that indeed only the first term in the delta 
function contributes. Neglecting all higher order terms in k/Pl, we obtain the amplitude 
squared for Fig. [2] as 

where {N^ — 1)/4N^ is the color-factor, and the hadronic Mandelstam variables s, t, and u 
are defined as s = {Pa + PbY, t = {Pa — PiY, u = {Pb — -Pi)^- In the limit Xb = 1, these 
variables can be related to the partonic Mandelstam variables as: s ~ s/xa, t ^ t/xa-, u ^ u 
at 0{Pj^, and we will have s + i + u = 0. Substituting the above results into the differential 
cross section, we obtain. 



c 



dyidy2dPld^q± AN^ 



9n m \ K^aMp + mqY + q]_ 



167r3 A^(x.,gi) 



K^h - 1) 



S2 t2 

(18) 



Comparing this result with Eq. ([71), and identifying the TMD parton distribution Xaf{xa, q±) 
defined in Eq. ([9]), we extract the partonic hard factor for the differential cross section from 
the diagram in Fig. 3 as 

_^Nl-12{l + v^ 

^qq'^qq' p a ATI 2o ' y^^) 



4Ar2 £2 



where s, t, u are the partonic Mandelstam variables defined before. This hard factor is just 
the well-known differential cross section da/dt for the partonic process qq' qq' . The hard 
factors for all other partonic channels can be calculated similarly, and they are all identical 
to the partonic differential cross section da/dt, and collected in Appendix A. 

Instead of calculating the differential cross section directly, we now show how to use the 



power counting techniques 3l|] to analyze the contribution from all partonic channels and to 
demonstrate the factorization of the hard factors. We then apply this approach to calculate 
the contributions to the spin dependent differential cross section. 
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FIG. 4: Factorization of the unpolarized cross section of Fig. into the TMD quark distribution of 
the proton and the hard factor represented by the partonic scattering Feynman diagram. 

We again first consider the simple example of the qq' qq' channel. The jet trans- 
verse momentum is the only hard scale, while all other momentum scales, such as the 
imbalance q± of the two jets and all other mass scales Mp, and A^, are relatively soft 
and much smaller than the hard scale. We expand the scattering amplitude squared in the 
small parameter q±/P±. Using power counting techniques, we identify the leading power 
contributions and match them to the momentum regions of the unobserved (or integrated) 
parton of momentum k'. In our current model, where a parton of momentum Pb scatters 
off a nucleon of momentum Pa, the leading power contribution is from the region when k' 
is parallel to P4. If k' is parallel to Pb, the contribution is power suppressed, because the 



denominator of the quark propagator (P4 — k') — m will be of order s = {Pa + Pb) ~ P 



instead of as is the case when k' is parallel to Pa- 

When the momentum k' is parallel to Pa, the quark propagator in Fig. [2] has a virtu- 
ality A'^{xa,qj_) ~ q'j_, which is much smaller than the hard scale Pj. Compared to the 
short distance physics in the partonic scattering taking place at the time scale 1/P±, the 
incident quark of momentum Pa — k' is long-lived, and the contribution from this subpro- 
cess can be separated into two parts: the parton distribution part relevant at the soft scale 
~ and the hard partonic part at the scale ~ 0{P±), representing the long distance 

and short distance physics, respectively. We demonstrate this factorization in Fig. IH where 
the contribution to the differential cross section from the process in Fig. [2] is factorized into 
a parton distribution multiplied by a hard factor. We can further separate the hard fac- 
tor into a partonic scattering amplitude squared stripped of any color, and the color-factor 
C„ = Tr(T'*T^)Tr(T"T'')/A/'2 = {N^-1)/AN^. In this way we can take into account the color 
decomposition of the subprocess in the case of unpolarized scattering. This is trivial for this 
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particularly simple partonic channel, but becomes a convenience for the more complicated 
ones. After subtraction of the parton distribution from the differential cross section, the 
hard factor in Fig. H] can then be written as 

Hqq'^qq' = "^^'^ X h^gi^gq' , (20) 

where a^ir/s^ represents a common factor and is introduced to simplify the notation, and 
hqqi^qq' = 2(s^ + u^)/i^. Thls of coursc agrccs with the above result in Eq. ( |T9l) . Most 
importantly, it demonstrates that the hard factor is model independent. It only depends 
on the partonic scattering cross section and the color-factor associated with the partonic 
diagram. It does not depend on the model we used to describe how the nucleon couples to 
the quark and diquark. 



B. Hard factors for the spin-dependent cross section 

A non- vanishing single transverse-spin asymmetry requires initial/final state interactions 
generating a phase. Because all initial and final partonic states relevant for dijet production 
are colored, both initial and final state interactions have to be taken into account for a 
complete result. In Fig. El we show all possible diagrams with gluon exchange between 
the diquark and the hard scattering part which may contribute to the SSA for the dijet- 
correlation, again for the case of an underlying qq' — > qq' process. The initial and final state 
interaction diagrams in Fig. [5]^a-c) contribute to a soft gluon pole, which corresponds to 
the quark Sivers function contribution we are considering in this paper. However, the gluon 
interaction with the internal gluon propagator shown in Fig. [S](d) does not contribute to a 



soft gluon pole, but to a soft fermion pole 



251]. This soft fermion pole contribution is not 



related to the quark Sivers function and will not be discussed in this paper. Therefore, in 
the following calculations, we will only consider the contributions from the diagrams (a)-(c). 

We will follow the same factorization approach that was used in last subsection for cal- 
culating the hard factors for the unpolarized cross section. We use the power counting 
method to factorize the additional gluon interaction between the hard scattering and the 
diquark into the Sivers parton distribution for every partonic channel. As discussed above, 
the leading power contribution is from the region of phase space where the outgoing diquark 
momentum k' is parallel to the polarized proton Pa- k'^ ~ P^, k'~ ~ (fll Pa ^^"^ ^1 ~ 
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FIG. 5: Initial and final state interactions between the scattering partons and the diquark of the 
polarized proton: (a) for initial state interaction, and (b,c) for the final state interactions diagrams, 
and (d) for interaction with the internal gluon propagator which does not contribute to the soft gluon 
pole corresponding to the quark Sivers function contribution to the dijet- correlation SSA. 

To obtain a leading power contribution, the gluon momentum / also needs to be parallel to 
Pa- This is so because the diquark propagator in Fig. [5] reads 

1 1 



(21) 



{k' + l)'^ + ie 2k'+l- + Oiql) ' 

which will be power suppressed unless l~ ~ (1±/Pa- This result can also be derived from the 
existence of a "pinch" singularity in the /~ integral over the quark propagator of momentum 
Pa — k' — I and the diquark propagator of momentum / + k' 35|. Since the momentum 
I of the gluon is parallel to Pa, its polarization will also be proportional to Pa, because 
the gauge-boson-scalar vertex is proportional to {2k' — Z)'^, as shown in Fig. [H Therefore, 
the gluon attaching to the initial- or final-state quark is longitudinally polarized, i.e., its 
polarization is along its momentum. Because of this, we can further decouple the gluon 
interaction with the external particles by using the eikonal approximation. For example, 
for the diagram with final-state interaction in Fig. Mjo), the gluon interaction part can be 
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FIG. 6: Eikonal approximation for the attachment of a longitudinally polarized gluon to external 
lines: (1-3) for the final state interactions when the gluon attaches to outgoing particles; (4-6) for 
the initial state interactions for incoming particles. The color-factor is not shown in the above and 
should read: T^- for (1,2,4,5) where a is the color index for the gluon with momentum I and i,j 
are indices for the quark lines; {—ifabc) for (3,6) where b,c are indices for incoming and outgoing 
gluons. 



reduced to 



u{P,)T%-zg)r 



(Pi - ly -ml + ie 



-2Pf ■ /+ + ie 
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u{Pi)T'' 



uiPAT" . . . 



(22) 



-/+ + le 

where the Gamma matrix 7" appears because of the interaction with a longitudinally po- 
larized gluon, and where a is the color-index for this gluon. In the above derivations, we 
have only kept the leading power contributions. 

Similarly, the initial state interaction in Fig. [5](a) and the other final state interaction in 
Fig. [5](c) can be simplified following the same eikonal approximation. The eikonal approx- 
imation will give rise to eikonal propagators and vertices. In Fig. [6l we summarize these 
eikonalized expressions for all the initial and final state interactions, including the attach- 
ments of the longitudinally polarized gluon to a quark or anti-quark line, or to a gluon line as 
well. The eikonal approximation for quark- and antiquark-gluon interactions with incoming 
or outgoing momentum can be derived following the above example. For the three-gluon 
interactions, we can derive the eikonal vertex and propagator by choosing physical polar- 
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izations for the external gluons. With this choice, we find the vertex and gluon propagator 
can be reduced to, for example for Fig. [6]^6): 

i-tg)Hfabc) [{I - PB)u'gp^ + (2Pb + l)pgt.u' + (-2/ - PB),gu'p] ^^^i ■ 

[-TB + I-) + 

~ Jp///)2+^e^-'^-^-) = _i+_^^9A-^fabc) , (23) 

where we have used the fact that the incoming gluon is transversely polarized, so that the 
polarization tensor can be chosen as 

pM pM' _|_ pp' pM 

E^A(^B)ef (Pb) = + ''bPa+PbPa (24) 

X Pa ■ Pb 

Similarly, we can simplify the longitudinally polarized gluon attachment to the outgoing 
gluon line as shown in Fig. [6]^3). 

These eikonal propagators and the associated vertices can be absorbed into the parton 
distribution for the polarized proton, after factorizing the diagram. The only difference be- 
tween this factorization and the one for the unpolarized scattering is the color-flow, because 
of the additional gluon exchange between the diquark and the hard part. The color- flow is 
sensitive to the difference between initial- and flnal-state interactions and depends on the 
speciflc partonic subprocess. Therefore, the generalized factorization would immediately 
fail if the eikonal propagators and associated vertices, which are absorbed into the TMD 
distributions, were sensitive to the process dependent color flows. However, at the leading 
order, one can always attribute the color-flow to the partonic hard part and keep the corre- 
sponding TMD distribution insensitive to it. Therefore, as compared to the unpolarized case 
(see Fig. Hj), for the single-spin-dependent cross section we shall obtain a more complicated 
factorization. For the case of the three diagrams of Fig. [5] the factorization is shown in 
Fig. [71 The contribution from each diagram can be factorized into the Sivers function and 
the hard partonic scattering amplitude squared, multiplied by a color-factor that takes into 
account the color-flow of the initial/flnal state interaction. However, in order to verify the 
generalized factorization beyond the leading order, one would have to show that additional 
initial-state as well as flnal-state gluon interactions could also be factorized into a sum of 
terms, each of which is factorized into a product of higher order TMD distributions and cor- 
responding color factors and partonic hard parts. For example, to verify the factorization 
at the next-to-leading order with two eikonal gluon interactions, one would have to show 
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FIG. 7: Factorization of the initial/final state interaction contributions to the SSA in dijet corre- 
lations into the quark Sivers function ( defined as in SIDIS), multiplied in each case by an adjusted 
color-factor and a squared scattering amplitude for the partonic process without color factor, similar 
to Fig. 

that the initial- and final-state interactions can be factorized into the same color factors 
and the lowest order hard parts as in Fig. 8, multiplied by the corresponding second order 
TMD distribution. We will leave the detailed study of the generalized factorization at higher 
orders to a future publication. 

At least to the leading order, the extracted color factors Ci with i = I , Fi, F2 in Fig. [7] are 
insensitive to the long-distance physics. That is, they do not depend on how the quarks and 
the gluon couple to the nucleon. This is because the factorization of the color-fiow between 
the hard partonic part and the TMD part results in a color configuration connecting an 
incoming and an outgoing quark (with color indices ij) and a gluon (with color a). The 
only possible color structure to describe this connection is T^", due to SU{3) invariance and 
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the Wigner-Eckart theorem. Therefore, one can reduce all the color matrices from the hard 
part to a simple color matrix T^j, and the coefficient will be exactly the color-factor in the 
factorization formula. In this way, the extracted color-factor cannot depend on the bottom 
part of the diagram describing the coupling of the quarks and the gluon to the nucleon, but 
will only depend on the color-flow in the hard partonic diagrams with an additional gluon 
insertion. For example, the color factors appearing in Fig. [7] can be derived in the following 
way, 

m ■■ ^MT'm{T''T% = ^7-. ^C, = ^y (25) 

CKb) : ^TT{T'mT'T'^T% = ^7™: ^ Cfi = , (26) 

1 /V^ _ 2 AA^ — 2 

m : -T:t{T'^T'^T%T'T% = -^n, ^ = , (27) 

where the color-matrices on the left side of the equations come from the partonic diagrams 
with one gluon attachment to the initial- or final-state quark line. The extracted color-factor 
Cj is for the initial-state interaction, and Cpi and Cf2 are for the final-state interactions on 
the lines with the momentum Pi and P2, respectively. 

In the factorization shown in Fig. [TJ the leading order Sivers functions for all three 
diagrams are the same. This is because we normalize all eikonal vertices to the one used 
in the calculations of the Sivers function in SIDIS (shown in Fig. [3]). The difference from 
the initial or final state interaction effects is summarized into the relevant color-factors. 
For other partonic channels, it is in some cases necessary to introduce an extra sign when 
the eikonal propagator contributes with an opposite phase compared to that in the Sivers 
function in SIDIS. For example, the eikonal propagators in the diagrams (2,5,6) of Fig. [6] 
contribute an opposite sign, whereas those of (1,3,4) contribute the same sign as in SIDIS. 
Because all the initial and final state interactions in Fig. [7] contribute the same sign, there 
is no sign change for their color-factors. 

Another important point is that the partonic scattering amplitude squared is the same 
for all the three diagrams in Fig. [71 It can be calculated directly from the relevant Feynman 
rules for the partonic diagrams and is also the same as the one for unpolarized scattering 
given above. After factorizing out the Sivers function from the cross section contributions 
for the diagrams in Fig. U\ the color factors for the three diagrams can then be summed to 
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give the final result for the hard factor for the qq' qq' subprocess: 



J^qq'^qq' " ^ + Fl + <-F2j tlqq'^qq' " J2 



(28) 



where hqq'^qqi has been defined above after Eq. (120|) . 

The above derivations can be extended to all other partonic channels involving the quark 
Sivers function. They can also be extended to the case of sub-processes initiated by the 
gluon Sivers function for which, however, more color-structures (like fabc and dabc for the 
three gluon coupling to the nucleon state) will emerge ^S^j- In this paper, we will focus on the 
processes with a quark Sivers function, which include all quark-quark and quark-antiquark 
scattering processes, such as qq' — > qq' calculated above, qq' gg', qq q'q' , qq — ^ qq and 
qq qq, as well as the quark-gluon process qg qg and quark-antiquark annihilation into 
a gluon pair, qq ^ gg. Following the above derivations for the qq' —>■ qq' process, the hard 
factors for the unpolarized and single-spin-dependent cross sections for these channels can 
be summarized by the following master formulas: 



where i labels an individual Feynman diagram (meant as the square of an amplitude or an 
interference between two amplitudes, see below) with the associated expression for the 
hard factor. is the color-factor for the unpolarized cross section, C} for the initial state 
interaction for the single-spin dependent cross section, and Cpi and ^r the final state 
interactions when the gluon attachment is to the lines of momentum Pi and P2 respectively. 

From the above analysis, we can easily see that the hard factors in Eqs. (1291) . pop are 
model-independent. The model-dependence is removed once we factorize out the unpolar- 
ized quark distribution or the Sivers function from the differential cross sections. From 
the above master formulas, the hard factors depend on the partonic scattering amplitude 
squared, and the color-factors, derived from partonic diagrams with or without the addi- 
tional gluon attachment. Both the color and partonic factors, and therefore the hard factors, 
are independent of the model for the proton used in our calculation. 

For the quark- (anti) quark scattering channels, we plot all the relevant diagrams in Fig. [8l 
The color-factors for unpolarized scattering can be straightforwardly evaluated. In the spin- 
depedent case, the color- factors for each partonic channel can be calculated following the 




(29) 




(30) 
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FIG. 8: Quark- (anti) quark scattering diagrams. The mirror diagrams of (3) and (6) are not shown; 
their contributions in Table I are identical. 



TABLE I: The color-factors and hard cross sections of individual diagrams of Fig. [8] for the quark- 
quark scattering channels. 
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above example for qq' qq', by reducing the strings of color-matrices for the partonic 
diagrams with an additional initial- or final-state gluon attachment to the simple T^j, and 
extracting the coefficients as the relevant color-factors. We list all these color-factors in Table 
I, including those for the unpolarized scattering. The color-factors in the spin- dependent 
case can also be calculated by contracting the string of color-matrices for a given partonic 
diagram with a T", analogous to what was done in 25|. These two methods yield the same 
answer. The partonic amplitude squared for each diagram in Fig. [8] can be easily calculated 
using the Feynman rules in Feynman gauge. We also list the results for the /i* in Table I. 
Combining the results given in Table I, we can obtain the hard factor for any particular 
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FIG. 9: Quark-gluon scattering diagrams. The mirror diagrams of (3), (5) and (6) are not shown; 
their contributions in Table II are identical. 

TABLE II: The color- and hard factors for the qg — > qg scattering channels in Fig. [9l where 
Cf = {N^c - l)/2A^c. 
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quark- (anti) quark scattering channel. For example, for the qq — *• qq subprocess, we have 
contributions from diagrams (1-3) of Fig. [HI So, the hard factor for this process will be 



Sivcrs 
qq^qq 



(1) 



.(1) 

'F2 



(2) 



.(2) 



cf ^ + c 



(3) 
Fl 



c 



(3) 
F2 



(31) 



where a factor 2 in the third term comes from the mirror diagram of (3). The full hard 
factors for all processes are listed in the Appendix. 

In Fig. [9l we show the diagrams for the quark-gluon scattering channel qg qg. Our 
results for this channel are listed in Table II. The color-factors are calculated following the 
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FIG. 10: The qq — > gg scattering diagrams. The mirrors of (3), (5) and (6) are not shown; their 
contributions in Table III are identical. 



TABLE III: The color- and hard factors for the qq — > gg scattering channels in Fig. (TOl 
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method described above. In order to factorize the cross section and to make the eikonal ap- 
proximation, as mentioned above, we choose a physical polarization for the external gluons, 
using the polarization projection in Eq. (12^ for the incoming gluon with momentum Pb, 
and similarly for an outgoing gluon, 

Y.'l{k)ei\k) = + ^''^^^+f'^^' , (32) 

A 

where k = Pi or k = P2. With these choices, we calculate the partonic amplitude squared, 
hgg^qg, for each diagram in Fig. [9] and list all results in Table II. Again, the full hard factors 
for this channel can be derived from our results in Table II using the master formulas in 
Eqs. ( 129|) and (!30|) for the unpolarized and polarized scatterings, respectively. 
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In Fig. [TOl we show the relevant diagrams for the qq gg channeL Our results for the 
color-factors and the squared partonic amplitudes are listed in Table III. The hard factors 
can be calculated accordingly, and their results are listed in the Appendix. 

We note that our approach for deriving the hard factors for the dijet-correlation can 
be straightforwardly extended to the case of jet-photon correlations in hadronic reactions. 



This process was also considered in Ref . 37|] . Here we have contributions from the partonic 
channels qg qj and qq —>■ g'J- The hard factors for both channels can be derived from the 
above results by replacing one final-state gluon by a photon. We list the full corresponding 
hard factors for both the unpolarized and the single-spin dependent cross sections, which 
agree with the ones of 37|], in the Appendix. It is interesting to note that the initial 
state interaction dominates for the qg — > qj channel, whereas for qq —>■ g'-f the final state 
interaction dominates. 

As the above calculations have demonstrated, at the order we have considered the TMD 
parton distributions extracted from the SIDIS (or Drell-Yan) process can be used to predict 
the cross sections or the azimuthal asymmetries for the dijet-correlation at hadron colliders. 
Of course, the initial and final state interactions contributing to the SSA in dijet-correlations 
will introduce new observable effects, reflected by the modified hard factors. For example, 
for the dominant channel qg qg for the dijet-correlation, the hard factor for the single- 
transverse-spin dependent cross section is about a factor 1/2 smaller than that for the 
spin-independent cross section, and the sign of this contribution is the same as that for the 
Sivers asymmetry in SIDIS. Such effects can be tested experimentally in the near future. 
Some of the phenomenological consequences of this have already been investigated in 38 1. 

In our analysis, we have adopted the SIDIS definition of the TMD parton distributions. 
We could also have chosen the definition according to the Drell-Yan process. The only 
difference between these two definitions is the direction of the gauge link. As is well- 
, the so-called naively time-reversal-odd parton distributions change sign 



known 
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between these two processes, whereas the time-reversal-even ones remain the same. So, if 
we change the definition of the TMD parton distributions to that valid for the Drell-Yan 
process, the hard factors for the single transverse-spin dependent cross section for the dijet- 
correlations will change sign as well. For the spin independent cross section, the hard factors 
remain the same. We stress that the physical (hadronic) cross sections do not change with 
the definition of the parton distributions. 
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We emphasize again that the hard factors that we have calculated in this section are only 
leading order. The small transverse momentum scale q±_ has been generated by the model 
TMD distributions that we have used. We have been able to factor out the partonic hard 
parts at (9(Pj_), keeping all g_|_-dependence in the TMD distributions. In the next section, 
we will take a different approach and treat the cross sections in collinear factorization, 
generating g_L perturbatively from (nearly collinear) gluon radiation. In this way we can 
to some degree investigate the role of perturbative QCD corrections, and their influence on 
factorization. We will find that the generalized factorization formulas in Eqs. ([SD,© remain 
valid at this order. 



III. COLLINEAR FACTORIZATION APPROACH 



In this section, we will use the collinear factorization approach to calculate contributions 
to both the spin-averaged and the single transverse-spin dependent dijet cross sections from 
partonic processes with one-gluon radiation. We will work in a kinematic region of "inter- 
mediate" imbalance between the two jets: Aqcd ^ ^ P±- When q± ^ Aqcd, we expect 
perturbative QCD and the collinear factorization approach to be valid for calculating both 
the q±- and the Pi-dependence. On the other hand, the limit q± <^ P± allows us to verify 
the generalized factorization formulas in Eqs. (I2]),(l5]) in the transition region where the TMD 
factorization could also be valid 28|]. By an explicit calculation, we will demonstrate that 
the contribution from one gluon radiation nearly parallel to the incident hadrons can be 
factorized into the same hard factors calculated in last section and the 0{g'^) perturbatively 
generated TMD parton distributions defined in Eq. (jl]), which effectively shows the validity 
of the generalized TMD factorization formulas in Eqs. ([2]) and ([5]) at this order and in this 
intermediate region. 

In Refs. 28| we performed similar calculations for the much simpler cases of the Drell-Yan 
and SIDIS cross sections, with the same result. This established that the two mechanisms 
for generating SSAs, one based on the Sivers function and TMD factorization, the other on 
twist-3 quark-gluon correlations and collinear factorization, are related and can in a sense 
be regarded as "unified". It is important to point out that the consistency between the 
TMD and collinear approaches in the intermediate region Aqcd ^ Q'± ^ -Pi is not sufficient 
to actually prove the TMD factorization formalism in the region q± ~ Aqcd- But, the 
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consistency is certainly a necessary condition that TMD factorization formulas like those in 
Eqs. ([2]) and ([5]) need to satisfy. It is the purpose of this section and this paper to show this 
consistency for the dijet momentum imbalance at the first non-trivial order in perturbative 
QCD. 

We will take the simple qq' qq' subprocess as an example to show how this works. 
The extension to all other channels will follow. For this channel the hard factors H^gi_,ggi 
and H^gJ^gg, for the unpolarized and polarized cross sections have been calculated in the 
last section. When the dijet-imbalance is large compared to Aqcd, "we can expand the 
factorization formulas, and the q± dependence will come from the TMD parton distributions 
at large transverse momentum. For example, Eq. ([S]) for the unpolarized cross section will 
become, 

dyidy2dPl<Pq^ Iaqcd«5x«Px 

= {xaq{xa, q±)xbq{xb) + Xaqixa)xbq{xb, g±) H } , (33) 

where Hg^^gg, is the hard factor, and where we only keep the two contributions from the 
incident quark distributions at large transverse momentum. q{x,q±) is the TMD quark 
distribution and q{x) the quark distribution integrated over transverse momentum. There 
could in general also be contributions from the soft factor and/or from fragmentation func- 
tions in the final state. In this paper, we are only interested in verifying the part of the 
factorization formula that involves the TMD parton distributions in Eq. (j4]). 

Because q± ^ Aqcd, the transverse- momentum-dependence of the TMD parton distribu- 
tions can be calculated in perturbative QCD, and the results will depend on the integrated 
parton distributions. This dependence is well known, see 3] for example, 

i^J_C f — 
27r2 (f^ ^ J X 



q{xa,q±) = tt^^Cf I —q{x) 



(34) 



where Cp = {N^ — l)/2A'"c, and q{x) is the integrated quark distribution as mentioned above 
and ^ = Xa/x. Note that we only keep the collinear part in the distribution^ i.e., 1 — ^ 7^ 0. 
It will certainly be necessary to consider the soft contribution at ^ = 1 [28]. Similarly, we 
can calculate q{xb, q±) at large q±. Substituting these two results into Eq. fl33l) . we obtain 

rfV"" , 1 a., 
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where ^ is as defined above and ^' = Xb/x' . In the next subsection, we will demonstrate 
by explicit calculation that the above result is reproduced by the direct calculation in the 
collinear factorization approach. 

Similarly, the Sivers function qT{xa,k±) at large k± can be calculated in perturbative 
QCD, and the result depends on the twist-3 quark-gluon correlation function Tp (the Qiu- 



Sterman matrix element) 25|, |40|: 
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where C is the proper gauge link to make the matrix element gauge invariant, and where 
the sums over color and spin indices are implicit. Keeping again only the contribution from 
the collinear part, we will have [28 ]. 
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(1-0 J - ' ^'(1-0- 

where as above ^ = Xa/x and Xg = {1 — ^)x. The definition of the quark Sivers function in 
Eq. (jlj) has been used to obtain the above result. We emphasize again that in this definition 
the gauge link is simple in the Feynman gauge, and goes to +00. As mentioned above, 
we have to take into account the perturbative expansion of the gauge link up to second 
order {0{g'^), see also the diagrams drawn in 128]). Substituting the above result and the 
unpolarized quark distribution q{xb,k±) of Eq. into the factorization formula Eq. ([2]), 
we obtain the single-spin dependent cross section in the intermediate transverse momentum 
region as 
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FIG. 11: Generic diagrams for quark-quark scattering contributions to the unpolarized and the 
single transverse-spin dependent cross sections in the coUinear factorization approach. 

In the above result, the term in square brackets comes from the Sivers function in the 
polarized nucleon at large- /c_l, while the last term corresponds to the large- /c_l unpolarized 
quark distribution from the unpolarized nucleon. We have also used the leading order 
relation llJl 



to derive the second term. 

In the following, we will derive the above spin-averaged and spin- dependent differential 
cross sections in the collinear factorization framework in the regime P± ^ q± ^ ^qcd, 
and verify the expressions in Eqs. fl35|) . fl38l) . This will provide a clear demonstration of 
the consistency between the collinear factorization approach and the TMD factorization 
formalism in the region of intermediate dijet imbalance. In Fig. [TTl we draw the generic 
Feynman diagrams for the unpolarized and single-transverse-spin dependent cross sections 
for the dijet-correlation in the collinear factorization approach. In these diagrams. Pa and Pb 
are the momenta of the two incident hadrons. For the single-transverse polarized scattering, 
Pa labels the polarized hadron. Pi and P2 are the momenta of the two jets in the final states. 
ka and kf, are the momenta of the two incoming partons, and for polarized scattering ki and 
/c2 will be needed to define the complete kinematics for the momentum flow associated with 
the polarized hadron. The blob in the center represents tree-level Feynman diagrams with 
the given initial- and final-state partons. In the collinear approach, jets are produced back- 
to-back by the Born diagrams, so that there will be no imbalance between the two jets if there 
is no additional gluon radiation. In order to obtain an imbalance, we have to add at least one 




(39) 
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gluon (momentum k') into the final state, as shown in Fig. [TTl The transverse momentum 
of k' will be equal to the imbalance of the two jets, i.e., k'j_ = —q± = —{Pi± + P2±)- 

We will use the simple qq' qq'g subprocess as an example to demonstrate the method 
and details of the derivation. All other partonic channels can be analyzed in a similar 
manner. In subsection A, we will perform the calculation for unpolarized cross section, and 
in subsection B, we will deal with the single transverse-spin dependent cross section. 



A. Unpolarized qq' — > qq'g 



In terms of collinear QCD factorization for unpolarized hadronic collisions [3l|], the con- 
tribution by the qq' qq'g subprocess to the dijet cross section can be written as (see 
Fig.ETa)): 

da^'^n r (jr dr' 1 

- ^ T^i^)^\^')T7r^^M,,^^^^^^^ (40) 



dyidy2dPl(Pqi_ J x x' ' ^ ^ 'IQs{2t^Y 

where s = {Pa + Pb)"^, l-Mqq'^qq'gl"^ is the spin and color averaged scattering amplitude 
squared for the partonic subprocess qq' qq'g, and q{x) and q'{x') are quark distributions 
for the incoming hadrons, A and B, respectively, at momentum fractions x and x'. In the 
collinear factorization approach, initial-state quark momenta are approximated as k^ = xP^ 
and fc^ = x'Pj^. Neglecting the invariant mass of the jet which is much smaller than the 
jet's energy, and keeping only leading powers in q±/P±, we can write the jet momenta as 
Pi = P± (^i, 1) and P2 = P± (^, -1). In Eq. (gO]) and the rest of this paper, the 
dependence on factorization and renormalization scales is suppressed. 

The basic qq' qq' process only has one Feynman diagram. Radiation of an additional 
gluon leads to a drastic increase in the number of diagrams. In Fig. [121 we show sample 
diagrams that are most relevant for our analysis in the limit P± ^ q± ^ Aqcd- We work 
in Feynman gauge, and classify these diagrams into different groups by using the power 



counting technique for this limit 3l|] . For example, if the momentum k' of the radiated gluon 
is nearly parallel to P4, only diagrams (1-5) contribute at leading power in q±/P±, and all 
other diagrams are suppressed by q±/P±. The contributions by these five diagrams will be 
factorized into the TMD quark distribution of the hadron with momentum Pa- Similarly, if 
k' is nearly parallel to Pb, only (5-9) contribute at leading power, and their contributions can 
be factorized into the TMD quark distribution of the hadron with momentum Pb- Diagrams 
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FIG. 12: Some of the Feynman diagrams for the qq' qq' channel contributing to the dijet 
imbalance in the coUinear factorization approach. Diagrams (1-5) contribute to the cross section 
at leading power when k' is parallel to Pa, which can be attributed to the parton distribution in 
hadron A, whereas (5-9) contribute when k' is parallel to Pb and belong to the parton distribution 
in hadron B. All other diagrams are power suppressed when k' is parallel to either Pa or Pb and 
will contribute to other factors in the factorization formula. 

(10,11) and other similar diagrams will contribute at leading power when k' is parallel to 
Pi, the momentum of one of the outgoing partons. Such contributions will become part of 
the jet and be subject to the jet definition. If instead of the jet a hadron is observed, these 
contributions belong to a TMD fragmentation function. Diagram (12) does not contribute 
to any of the above cases, but will contribute when k' becomes soft and thus give rise to 
part of the soft factor in the factorization formula. As mentioned above, in this paper we 
will study the factorization of the coUinear gluon interactions into the gauge invariant TMD 
parton distributions. We will therefore focus on the contributions when k' is parallel to either 
Pa or Pb. The other contributions will be important to fully understand the factorization 
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formalism, but we leave them to a future publication. 

In the following, we will work out the result for diagrams (1-5) in some detail, in order to 
see how their contributions can be factorized into the TMD quark distribution of hadron A. 
We choose the kinematics so that the momenta of two incident hadrons only have light-cone 
"plus" and "minus" components, respectively: Pa = (P4,0^,0_l) and Pb = (0+,Pg,0_L). 
Then, in the collinear approximation, the momentum ka of the initial-state parton a only 
has a plus component while has only a minus component. For extracting the leading 
contribution in q±/P±, it is convenient to decompose the total dijet momentum q = Pi + P2 
and the momentum k' of the "unobserved" final-state gluon into 

q = Cka + C' kb + q± , 

k' = {i-^)ka+il-ah + k'^ , (41) 

— * 

where k'j_ = —q±, and ^ and ^' are the fractions of the initial quark momenta carried into the 
dijet final state. The partonic Mandelstam variables are given at leading power in q±/P± by 

s = i^ka + e'fcb)' = a'2fc, ■ k, = 2Pi ■ P2 , 

i = i^k, - Pi)2 = -^2ka ■ Pi = -e'2fc, ■ P2 , 

U = i^ka - P2)' = -^2ka ■ P2 = -^'2k, ■ Pi , (42) 

where ^ = Xa/x and ^' = x^/x' with Xa = ^ (e^^ +6^^), Xf, = ^ (e~^^ +e~^2), and s = XaX^s. 

When k' is parallel to Pa and k'j_ ^ Pi, the on-shell condition k''^ = (1 — ^)(1 — C){2ka ■ 
kb) — (^j_)^ = leads to (1 — ^') ^ while (1 — ^) remains large. The delta function in 
Eq. (140|) can then be reduced to 

^((m— f^Y^- (43) 

With this approximation, the squared amplitude will become much simpler. For example, 
the contribution by diagram (1) will be 



where the color-factor has not been included. For this diagram, it is very simple, and can 
be easily factorized as 

^Tr(T"T"T^T'=)Tr(T''T^) = Cf x ^Tr(T''T^) Tr (T*T") = x C„ , (45) 
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FIG. 13: Factorization of the contributions from diagrams (1-5) into the TMD quark distribution 
of hadron A times the hard factor, when the momentum k' of the radiated gluon is nearly parallel 
to Pa ■ These five diagrams are the leading ones among those in Fig. \1SX 



where C„ = {Nl — l)/ANl is the color-factor for the partonic process qq' — > qq' . Substituting 
the above results into Eq. pUj) . we obtain the contribution to the differential cross section 
from this diagram as 

2is^ + u' 



J uu 

dyidy2dPl<PqA 



iFiglMl) 



t2 



Xbq{xb) I dxq{x)^CF-^^{l - ^ (46) 



where the factor in the square brackets is the hard factor H^^,_^gg, for the TMD factorization 
formula, Eq. (jl]), and is the same as that calculated in the last section in Eq. f|T9|) . Other 
than Xbq{xb), the rest of the right-hand-side of the above equation is a part of the TMD 
quark distribution in Eq. ( IMl) . We illustrate this factorization in the upper panel of Fig. [T3l 
Diagrams (2-5) are more complicated, but their calculations are straightforward using 
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the above techniques. First, their color factors can be formulated as follows: 

(2) : ^Tr(T"T^r"T^)Tr(T^T'=) =CfxCu + -^'-fabcMT^^T'T') , 

(3) : -^Hfa,c)TTiT^T''r)TT{TV) = --ll/,,,Tr(r«T''T^) , 

(4) : -i^Tr(T"T^r")Tr(r'^T"T^) , 

c 

(5) : ^Tr(T"T^T'^)Tr(T'^T^T'=) , (47) 

where we notice that the only difference between the color factors for diagrams (4) and 
(5) is the T'^T^ vs. T^T'^ in the second trace. After applying Eqs. (H2|) . (143|) . the squared 
amplitudes (without the color-factors) for these diagrams are reduced to the following simple 
expressions: 

+ 1 



\M\h 


= 


\M\h 


= 


\M\U 


= 


\M\h 


= -1 



6 



6- 



l-M(4)P, (48) 

where the contributions from diagrams (4) and (5) are equal and opposite. Thus, when color 
is included, their sum will depend on the difference of the color-factors in Eq. (H7I) . 

l^Tr(T''T''T'm(T^rT^) - -^TT(T''T^T^m(T''T^T'') = — ^-fabcTriT^'T^r) . (49) 

2 

From the above results, we find that the contributions from diagrams (2-5) only have two 
color-factors: Cp x and —ifabc'~^T^{T"-T'^T'^)- However, it is straightforward to see that the 
contributions coming with the color structure —ifabc^'^{T°'T^T'^) cancel each other, so that 
the only contribution will have the color factor Cp x C„. Therefore, the final result for the 
amplitude squared from diagrams (2-5) will be, including the color-factor. 



1-^1(2-5) = <-i^ X p2 y^^) 



Inserting the above into Eq. (HOi) . we find the contributions by diagrams (2-5) to the differ- 
ential cross section: 
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where H^gi^ggi is the hard factor for the TMD factorization formula mentioned above. 
Clearly, apart from Xbq{xb), this result can also be interpreted as part of the TMD quark 
distribution for hadron A, multiplied by the hard factor. We display this factorization in the 
lower panel of Fig. [121 It is important to note that we have taken into account all possible 
attachments of the radiated gluon to the hard part to obtain the right answer. Especially, 
the diagram with the gluon attached to the gluon propagator is not power suppressed, and 
its contribution is crucial to achieve the final factorization. Summation over all gluon at- 
tachments leads to a gauge link contribution as shown in Fig. [131 Of course, this result is 
also consistent with the Ward identity that allows to factorize all gluon attachments into a 
single gauge link. 

After summing all the contributions by diagrams (1-5), including the mirror diagrams 
of (2-5), we obtain the final result for the differential cross section for the dijet-correlation 
when the radiated gluon momentum k' is nearly parallel to Pa'- 



dyidy2dPl(P'q^ 



2ti^ q]_ 1 — 4 



(52) 



= Hg^,_^gg,{s,i,u) [x b Q ^ dxq^x) 

k'ocPA 

which reproduces the first term in the TMD expansion result in Eq. f[5^ . 

Similarly, from diagrams (5-9) and the corresponding mirror diagrams, we evaluate the 
leading contribution in q±/P± when the gluon momentum k' is nearly parallel to Pb- Com- 
bining with Eq. ( !52l) . we find: 



^^(qq') 



dyidy2dP]_d'^q_i 



coll. fac. 



Q_\_ ^ J 



Px>gx>AQCD 



(53) 



:i-o ' (i-e 

with ^ = Xa/x and ^' = Xb/x' as introduced above. We notice that in Eq. f[S51) the physics 
at the scale q±_ is decoupled from that at Pi, reproducing the first order TMD factorization 
result in Eq. (1351) . 

In summary, by explicit calculations of the Feynman diagrams for the partonic process 
qq' — i> qq' g which contributes to the dijet-correlations in hadronic reactions, we have demon- 
strated that when the momentum of the radiated gluon k' is almost parallel to either P^ or 
Pbi the contribution to the differential cross section in Eq. fl53p can be factorized into the 
perturbatively generated TMD quark distributions defined in Eq. (jlD for the two incident 
hadrons, multiplied by the same hard factor at C(-Pl) as calculated in the last section using 
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the lowest-order 2-^2 Born diagrams. This can be viewed as support for the generahzed 
factorization formahsm in Eq. ([5]) for the dijet-correlations at hadron colhders. 

B. Single transverse-spin dependent cross section 

In this section we calculate the leading contribution to single transverse-spin dependent 



cross section Aa{S±) in the limit P± ^ q± ^ ^qcb in the twist-3 ETQS approach (24, l25 |. 
The difference between the unpolarized cross section calculated in the last subsection and 
the transverse-spin dependent one is that the latter involves an additional gluon from the po- 
larized proton, which interacts with partons in the hard part, as shown in Fig. [TlTb). When 
both Pi and q± are much larger than the typical hadronic scale Aqcd, the collinear factor- 
ization approach should be valid for describing the SSA [39], and a nonvanishing Acr(5'_L) 
is generated by the ETQS mechanism in terms of twist-three transverse-spin dependent 



quark-gluon correlation functions 2J, |25|]. In the ETQS formalism, the contribution of the 



subprocess {g)qq' — > qq'g to Aa{S±), shown in Fig. fTTTb). is generically given by 

c^r^yli^Pj^"^ = / ^ d-^d-^ TFixu X,) q'ix') ^^-i-^ SiikT) n,,,,,^,,, , (54) 

where 7i represents a partonic hard part, xi and X2 are the momentum fractions of the quarks 
from the polarized hadron A on the two sides of the cut shown in Fig.[TT](b). and Tp(xi, X2) is 

uon correlation function defined in Eq. fl36|) . extracted 



the corresponding twist-three quark-i 
from the lower blob in the figure [25| 



3- 

The strong interaction phase necessary for a non-vanishing Aa{S±) arises from the in- 
terference between the imaginary part of the partonic scattering amplitude with the extra 
polarized gluon {kg = XgPA) and the real scattering amplitude without a gluon in Fig. [TlTb). 
The imaginary part comes from taking the pole of parton propagator associated with the 
integration over the gluon momentum fraction Xg. For a process with two physical scales, 
Pj_ and g_L, tree scattering diagrams in Fig. [TlTb) have two types of poles, corresponding 
to a;^ = ("soft-pole") \2d\ and 7^ ("hard-pole") 28|. When calculating the partonic 
scattering amplitudes, we have to attach the polarized gluon to any propagator of the hard 
part contained in the light circles in the diagram of Fig. [TlT b). If the polarized gluon at- 
taches to the external quark lines either in the initial state or in the final state, the on-shell 
propagation of the quark line will generate a soft gluonic pole. A hard pole arises when 
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internal quark propagators go on-shell with nonzero Xg. In Figs. [T3] and [15] we show the 
diagrams for the soft pole contributions, and in Figs. [T6lfT8] the ones for the hard pole contri- 
butions. In these figures, we only show the diagrams with the additional gluon attaching to 
the left of the cut line. Their mirror diagrams for which the gluon attaches to the right are 
not shown, but are included in the final results. Certainly, because of the additional gluon 
attachment, we will have many more diagrams to calculate, compared to the calculations in 
the last subsection for the unpolarized cross section. However, again, by using the power 
counting technique, we are able to group all diagrams according to, for example, whether 
the radiated gluon's momentum k' is nearly parallel to the polarized hadron's momentum 
Pa, or to that of the unpolarized one, Pg. Similar to the last subsection, we will focus on 
these two types of contributions, because we want to investigate how they may be factorized 
into the perturbatively generated TMD parton distributions defined in Eqs. flMl) . fl37|) . All 
other contributions are either associated with final-state effects, or with the soft factor. In 
the evaluations of the scattering amplitudes, we always keep only the leading power contri- 
butions in q±/P± and neglect all corrections that are power-suppressed. In this way, we can 
clearly investigate how the factorization works, and what the appropriate definition of the 
TMD parton distributions is in this case. 

The calculations of the soft-pole and hard-pole contributions to the single spin asymmetry 
for the dijet-correlations follow the same procedure as we used for the SIDIS and Drell- 
Yan processes 28|], because they correspond to the same kinematic limit, the intermediate 



transverse momentum region. In Ref. [28j, we first calculated the full differential cross 
section valid for all transverse momenta with q± ^ Aqcd, including the region q± ^ Q 
(where q± is the transverse momentum and Q the virtuality of the photon in SIDIS or the 
Drell-Yan process). The results were then expanded in terms of q±/Q to obtain the leading 
power contribution in the intermediate transverse momentum region. As mentioned above, 
however, in this paper, we will utilize the power counting technique from the beginning 
in order to simplify our calculations. This means that our results are only valid in the 
intermediate transverse momentum region Aqcd q± P±- A full calculation of the 
qT dependence of the dijet cross section also at g_L ~ Q would be extremely tedious and 
not really provide any additional insights. As a cross check, we have also recalculated the 
Drell-Yan SSA in the intermediate transverse momentum region following the method used 
in this paper, and found results identical to our earlier ones in 281]. 
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FIG. 14: Soft pole contributions when k' is parallel to Pa- The mirror diagrams where the polarized 
gluon attaches to the right hand of the cut line are not shown but included in the final results. 

The method for calculating the single transverse-spin asymmetry for hard scattering pro- 
cesses in the twist-three approach has been well developed and documented in the literature 



25 



28 



40 



4l|. In the following, we will just outline the main steps of the calculations 



and highlight the unique features for this particular 



prob 



28|, 



4d, 



em. For further details of the 



4l| . The coUinear expansion is 



twist-three calculations, we refer the reader to [25|, 
the central step in obtaining the final results. We perform our calculations in a covariant 
gauge. The additional gluon from the polarized hadron is associated with a gauge potential 
A'^, and one of the leading contributions comes from its component A'^. Thus, the gluon 
will carry longitudinal polarization. The gluon's momentum is dominated by XqPa + kg±, 
where Xg is the longitudinal momentum fraction with respect to the polarized proton. The 
contribution to the single-transverse-spin asymmetry arises from terms linear in kg± in the 
expansion of the partonic scattering amplitudes. When combined with A~^, these linear 
terms will yield d-'-A'^, a part of the gauge field strength tensor F-^~^ in Eq. (136|) . In this 
coUinear expansion, we further keep kg± <^ q± because q± is a relative hard scale compared 
to kg±. Since kg = k2 — ki, the kg± expansion of the scattering amplitudes can be performed 
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FIG. 15: Soft pole contributions when k' is parallel to Pb- 



through the transverse momenta of ki and which we can parameterize in the following 

way, 

ki = xiPA + ki^, k2 = X2PA + k2±, (55) 

where we have neglected the minus components of the momenta, as these do not contribute 
to the linear terms in the expansion in kg±. From momentum conservation, we know that 
kg± = — ki_i_. Therefore, the collinear expansion in kg_i_ can be replaced by expansions 
in ki± and k2±: 



H 



dT~L dT~L 



where TC represents the amplitude squared for the partonic process {g)qq qqg including 
the delta function for the phase space integral of k', 6{{k')^) in fl54|) . Because of gauge 
invariance, the terms linear in ki± and k2± can always be combined after summation over 
all diagrams, 

(57) 



^p-|-|fclX=fc2±=0 



on 
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FIG. 16: Hard pole contributions when k' is parallel to Pa- 

The final expression for tlie collinear expansion can tlien be written as 

d'hi 

'H(g)qq^qqgikl±, k2±) = TL(^g)qq^qqg{Q , 0) + | fcix=fc2X=0 • (58) 

In order to obtain the complete result for the right-hand-side of the above equation we have 
to keep track of the momentum flow in the hard parts. One important contribution of the 
kgi_ expansion comes from the on-shell condition for the radiated gluon, whose momentum 
k' depends on kg^. This leads to a term involving the derivative of the correlation function 
Tp. In addition to the derivative contributions, also non-derivative terms can arise from the 
/cg_L-expansion of other parts of the scattering amplitudes. In the following, we will calculate 
the derivative and non-derivative terms separately. 



1. Derivative Terms 

The derivative contribution is simpler. As summarized in [25!], the derivative terms come 
from two parts of the collinear expansion. One is the on-shell condition for the radiated 
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FIG. 17: //arc? po/e contributions when k' is parallel to Pb- 

gluon, 5{{k'Y), mentioned above, and the other is the double pole contribution in final state 
interactions. In the derivation of these derivative contributions, we only need to focus on 
these two parts in the kg±_ expansion of the partonic amphtudes Ti and can neglect all other 
contributions. 

We first discuss the contributions by the on-shell condition {k'Y = 0. Since the kg± 
momentum fiow will be different when the polarized gluon attaches to the left or to the 
right of the cut, we denote the radiated gluon's momentum as fc^ or k'j^ in these two cases. 
From momentum conservation, we find in Fig. [TH for example, that k'j^ can be written as 

k'L = h + kh-Pi-P2 = X2PA + k2± + kb-Pi-P2^k' + k2i_ , (59) 

where k' = xPa + x'Pb — Pi — P2- Here, we have also approximated X2 by x; their differ- 
ence will contribute to the non- derivative terms but not to the derivative ones. With this 
decomposition, (A;^)^ becomes 

{k'^y ^ {k' f - 2k2± ■ (Pix + P2±) = {k'f - 2k2x ■ q± , (60) 

where we have used the relation q± = Pi± + P2±. Here we have also neglected higher-order 
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FIG. 18: Power suppressed hard pole contributions when k' is parallel to Pb- 

terms in A;2_l which do not contribute to the single-spin asymmetry. Similarly, when the 
gluon attaches to the right, we find for the on-shell condition: 



(61) 



It is easy to see that the above two expansions differ only by kgj_ ■ qj_. We know that the 
gluonic poles have opposite signs for the two attachments, and that apart from this difference 
the squared amplitudes for them are identical. So we can combine the two contributions 
and find the following result for the expansion of the delta function: 

d6{{k'y) 



5{{k',Y)~5{{k'^Y) = -2k,^-q^ 



(62) 



d{k'Y ' 

We can further rewrite the derivative of the delta function as a derivative with respect to x, 

dm?) ^ d5{{k'Y) 1 , . 

d{k'Y dx d{k'Y/dx ' ^ ' 

By using the momentum decomposition k' = xPa + x'Pb — Pi — P2, the derivative d{k'Y/dx2 
can be written as 



d{k 



.'\2 



dx 



2k' -P . 



(64) 



This leads to an important observation: the derivative contributions will be power suppressed 
when the radiated gluon's momentum k' is parallel to the unpolarized hadron's momentum 
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Pb- This is because if k' is parallel to Pb, '2k' ■ Pa will be of order s ^ P]_, and from Eq. (163|) 
we can easily see that the contribution will be suppressed by 1/P]_- However, if k' is parallel 
to Pa, 2k' ■ Pa will be of order q'j_, instead of order P^, so that in this case the derivative 
contribution is not power suppressed. We can then further parameterize k' as 

^'^-(^-o^p^ + ^,^J^p^,p n-,i, (65) 

where the coefficient in front of (the minus component of k') comes from the on-shell 
condition {k'Y = 0, and where we have also used the relation fcf = —q'^. Substituting the 
above into Eq. (!64|) . we find that 

dx x(l - ■ ^ ^ 

This expression will appear in the final result for the derivative terms. 

When the polarized gluon attaches to the final state particles, there will be double pole 
contributions to the single-spin asymmetry [25]. This is because the position of the pole itself 
will depend on kg±, and the expansion of the associated propagator will lead to a double 
pole. For example, in the final state interaction diagram of Fig. fT^ 6i ). the propagator giving 
a pole reads 

(67) 



(Pi - + ze (Pi - XgPAy - 2P,^ ■ V + ' 
The expansion in kgj_ w ill g ive a double pole. This double pole will make a contribution 



to the derivative terms j25j. However, it turns out that this double pole contribution to 
the derivative terms is power suppressed in the intermediate transverse momentum region 
AocD -C (?± ^ P±. In order to demonstrate this, we will follow the method described in 



40l | to calculate the double pole contributions to the derivative terms. We recall that it is 
the imaginary part of the above propagator that contributes to the single-spin asymmetry. 
We can first evaluate the imaginary part and perform the kg± expansion afterwards. The 
imaginary part is proportional to the pole of the propagator in fl67|) . 5((Pi — kgY), for which 
the kg± expansion gives 

6{{Pl - kgf) = 6{{Pl - XgPA?) + 5'((Pl - XgPAf){-2Pll. ' k g ^) . (68) 

For the derivative contribution, the double pole will have the same form as that for the 
expansion of the on-shell condition for the radiated gluon discussed above. When we convert 
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the derivative of the delta function in Eq. (l68ll to a derivative with respect to Xg, we find a 
suppression factor proportional to 

Therefore, this double pole contribution is power suppressed in the q±/P± expansion. Sim- 
ilarly, the double pole contribution to the derivative terms from the final state interaction 
on the line with momentum P2 like the one shown in Fig. [TD^ci) is also power suppressed. 

Another important observation is that the derivative term in the hard pole contributions 
vanishes when all diagrams are summed. This has been shown in the explicit calculations 



for the Drell-Yan and SIDIS processes 



281], and has also been demonstrated from a more 



general argument based on the analytic property o: 



invariance of the twist-three matrix elements 4l|]. Here, we will give another proof for 



the momentum expansion and on gauge 



this result based on the Ward identity. Take one of the hard pole contributions shown in 
Fig. [T6](a) as an example, where all three diagrams (ai-as) contribute to the same hard pole 
as indicated by the bars in the propagators. As discussed above, the derivative terms come 
from either the expansion of the delta function for the on-shell condition for the radiated 
gluon, or from the double pole in the final state interaction (which also corresponds to a 
delta function as we showed above). When we expand the delta functions, all other parts 
of the amplitudes will be evaluated at kg± = 0. Comparing these three diagrams ai-os, 
we find that the only difference exists in the vertices where the polarized gluon attaches to 
the incoming quark, outgoing gluon, or the quark propagator, while the rest of the squared 
amplitudes is the same in each case. Thus, we can separate off this vertex part, and sum 
the remaining parts of the diagrams. These three vertex parts form a gauge-invariant set of 
lowest-order qg qg scattering diagrams with all external partons on mass shell. Without 
kg±-fiow in these vertices, the attaching gluon only has longitudinal momentum {kg = XgPA), 
and its polarization is also along the longitudinal direction. From the Ward identity, the 
sum of these three diagrams vanishes because all other partons in this lowest-order scattering 
amplitude are on mass shell and gluons have physical polarization. The situation is described 
in Fig. [IS 

In summary, in the intermediate transverse momentum region Aqcd q± <C Pl, we 
will only have leading-power contributions to the derivative terms from the expansion of the 
delta function for the on-shell condition of the radiated gluon k', and when k' is parallel to 
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FIG. 19: There is no contribution to the derivative terms from the hard poles because of the Ward 
identity. 

TABLE IV: The color factors for the diagrams in Fig. [TH 



(1) 



(2) 



(3) 



(4) 



(5) 



1 1 

2Nc 



1 1 

2Nc 2M 



_2 3_ 



_2 i_ 

" 2Wc 2Nc 



1 1 
2Nc ilVj 



1 jVg+1 
2Nc^M' 



1 1 
'2Nc 4 



1 

'2Nc InT 



1 1 

2Nc2m 



1 N^-2 
'2Nc 1W 



1 N^-2 
'2Nc 4m 



]_ 

2Nc 2Nc 



2Nc 2Nc 



Pa- We have shown these diagrams in Fig. [TH This part of the contribution is relatively 
easy to derive, as we mentioned above, by multiplying a factor (1 — to the unpolarized 
cross sections. The only additional involvement is the color factor. We list the color-factors 
for all diagrams of Fig. [THin Table IV. Each of the color factors has a factor l/2Nc. This 
can be easily seen for the color factors for diagrams ai, bi and ci, for which we have the 
following simplifications, 

1 



ai 
bi 

Cl 



1 1 
1 1 

1 1 



Tr T' T"^T"^) Tr (T'"T"^T"^) Tr ^2~"'T"^T"^) 

Tr (T'''T"^T'''T"^T"^T"^)Tr(T"^T"^)Tr(T"^T"^) 

Tr ^2~'''T"^T'''T"^T"^) Tr ^2~'"T"^T"^) Tr ^2~"^X"^T"^) 



2Nr 



X Ci 



2Nr. 



X C 



2Nr. 



X C 



F2 



(70) 



where Cj, Cpi and Ci?2 have been defined above in Eqs. fl25|) . fl26l) . fl27j) . Combining these 
color factors with the squared amplitude calculated in the last subsection for the unpolarized 
cross section, we find that the contributions to the derivative terms can be factorized into 
the hard factor given in Eq. (1281) times the derivative terms of the perturbatively generated 
Sivers function for the polarized nucleon in Eq. fl57|) . We illustrate this factorization in 
Fig. [201 For example, the derivative contribution from diagram (oi) is factorized into a soft 
pole diagram contribution to the derivative term in the Sivers function where the radiated 
gluon attaches to the quark line shown in Fig. [201 multiphed by the color-factor C/ and 
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FIG. 20: Illustration of the factorization of the soft pole contributions when k' is parallel to Pa- This 
part can be factorized into a soft pole contribution to the perturbatively generated Sivers function 
in polarized hadron Pa- This factorization applies both to the derivative and the non- derivative 
contributions from the soft poles. The circle with the attached gluon in the Feynman diagrams 
represents all possible gluon attachments. Likewise, the gluon attachment to the circle in the Sivers 
function diagram includes the diagrams with gluon attachments to the gauge link and to the quark 
line. 



the partonic amplitude squared for qq' qq'. hi and Ci will be factorized into the same 
diagram contribution to the quark Sivers function, but with different color-factors Cpi and 
Cf2-, respectively. The same conclusion holds for all other diagrams. For example, the sum 
of diagrams (02-05) is factorized into a soft pole diagram contribution to the Sivers function 
where the radiated gluon attaches to the gauge link, multiphed by the color factor C/ and 
the partonic amplitude squared for qq' — > qq' - 

Finally, we present the total result for the derivative contribution by the qq' process to 
the single-spin asymmetry for the dijet-correlation, 

k'ocPA 



dyidy2dPld?qA 



a/3 Ca J3 



QQ -^qq 



as 
2^ 



1 



2Nr. 



x-^Tf{x,x) 



[1 + e] , (71) 
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where the superscript "D" indicates the derivative term, and where the single-scale partonic 
hard part H^^^^^^^, has been defined in Eq. f l28|) . The result in Eq. flTTl) indeed reproduces 
the derivative terms in the first order TMD factorization formula in Eq. (1381) . 

2. N on- derivative terms: k' parallel to Pa 

In addition to the derivative contributions, the soft poles also contribute to the non- 
derivative terms. One way of calculating the non- derivative terms is to follow the kg±_-?low 
in the scattering amplitude squared, summing up all contributions. Keeping in mind that 
kg± = k2± — ki±, we have to keep the full dependence on k2± and ki± in the calculations. It 
is easy to keep track of the explicit ki± dependence coming from momentum conservation. 
In addition, we also need to include the "indirect" kij_ dependence resulting from kinematic 
constrains: the on-shell condition of the radiated gluon and the double poles in the final 
state interactions. For example, the on-shell delta function for the momentum k' of the 
radiated gluon requires that the longitudinal momentum component depends on kg±, which 
will lead to a kg± dependence for the longitudinal momentum fractions of ki and k2 as well. 

As mentioned above, if the polarized gluon attaches to the left side of the cut, the radiated 
gluon will have momentum k'j^ = k2+ki,—Pi—P2- Since k2 has transverse momentum k2±, the 
transverse component of k'^ will be equal to k'ij_ = k2± — Pi± — P2± = k2± — q±. Furthermore, 
the minus component of A;^ is fixed by momentum conservation to k'j^ ^ A;^ — — , 
where the contribution by fc^ has been neglected, because it is of order 0{k1j_). Therefore, 
k'^ does not depend linearly on kg_i_. The plus component of fc^, on the other hand, will 
have linear kg± dependence because of the on-shell condition. We therefore parameterize A;^ 
as 

kt = (1 - Ox (1 + 0{k2±)) Pj^ + f3P^ + ki^^ - q^^ , (72) 

where we have indicated the term linear in k2i. in the plus component of fc^. As mentioned 
above, the minus component [3 does not depend linearly on kg^ and can be determined to 
be /3 = q]_/2{l — C,)xPa ■ Pb- By using the on-shell condition (fc^)^ = 0, we can then obtain 
the kg ^-dependent term in the plus component of /c^, and the final parameterization result 
for k'l^ is 

?). (l - ^) Pi - + - . (73) 
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It is easy to check that {k'^)'^ = up to terms quadratic in k2±. Prom this result, we can 
also determine the A;2±-dependence of the plus components of ki and k2- For example, k2 
can be directly calculated from momentum conservation, and we find 

2q± ■ k2±' 



^2 



X 



i-(i-0- 



P'i + k(i 



2± 1 



(74) 



where the minus component of k2 has again been neglected. We next obtain the parameter- 
ization for ki by ki — k2 — kg, where kg is determined by the pole and will be different for 
different poles. For the initial state interaction, the pole {Pb -\-kgY — gives k'^ — kgj_, and 
the longitudinal momentum fraction does not depend on the linear term in kg±. However, 
for the final state interaction with gluon attachment to the line with momentum Pi, the 
pole will be given by (Pi + kg)"^ = 0, and the plus component of kg will depend on kg± 



through fc^ 



Pi-k„ 



9 Pa-Pi " ^ ' '"ff-L 

P2-k. 



P^ + Similarly, the final state interaction with P2 will lead to 



p^^P^ + Zc^j^. In this way, the parameterizations for ki and kg will vary for different 
poles, whereas that for k2 will always remain the same. This is because the parameterization 



of k2 comes from the on-shell condition for the radiated gluon, which is not sensitive to the 
location of the poles. 

When the polarized gluon attaches to the right of the cut lines, on the other hand, it is 
ki that is determined by the on-shell condition for the radiated gluon, (A;^)^ = 0. Following 
the same method as above, we find 

2q± ■ ki±_ 



kt = (1 



ki ^ X 



0^1 



Q± • Q± 
2q± ■ ki± 



-n 



2(1 - OxPa ■ Pb 



^A ^ '^1± 



(75) 



Q± ■ q± 

Furthermore, in this case k2 — ki + kg, with kg determined from the position of the poles; 
k2 will be different for the different poles. 

With the above parameterizations for all the external particles' momenta, the calcula- 
tions for the linear expansion in kg±^ are straightforward. The final result for the soft pole 
contributions to the non-derivative terms when k' is nearly parallel to Pa have the following 
form: 

_ rrSivers ^ ^±y± 

k'(xPA 



dyidy2dPl_d?qA 



[xb q'{xb)] 



X 



Ols_ 

27r2 



2Nr 



I 



dx^TF{x, x) 



2e - 



(76) 
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where H^'^^l^, is given in Eq. (l28l) . The separate factorization for initial and final state 
interactions diagrams in Fig. [20], which we used to describe the factorization of the derivative 
contributions, also applies to the non-derivative ones. For example, the contributions from 
diagrams (ai-as) of Fig. [T3] can be factorized into the Sivers function shown in Fig. [2Ul 
multiplied by Cj and the partonic scattering function, including both the derivative and 
non-derivative contributions. The same happens for the final state interaction diagrams 
(61-65) and (C1-C5). 

There are also hard pole contributions when k' is parallel to P4. We have shown these 
diagrams in Fig. [161 As discussed above, the hard pole diagrams only contribute to the 
non-derivative terms. To calculate their contributions, we follow the above method as for 
the soft pole contributions. Because the hard poles appear in the internal propagator, we 
have to take all possible gluon attachments into account. For example, diagrams ai-a^ of 
Fig. [16] display all contributions with the same hard pole (indicated by the short bars in the 
propagators), with three different attachments of the polarized gluon. 

Another important feature for the hard pole contributions is that, after taking the pole, 
the polarized gluon will have nonzero longitudinal momentum fraction, i.e., 7^ as we 
mentioned above. This is in fact the reason for referring to these poles as "hard" poles. 
For example, the hard poles in the diagrams 01-03 have the momentum kb — k' + kg with 
kg = XgPA, and their pole arises at the following value for Xg-. 

-{h-k'f ^ 2h-k' 
2Pj,.{k,-k')^ 2k,.PA~^ ^^""^ ^ ^ 

where the last equality comes from the fact that k' is parallel to Pa and is dominated by its 

plus component as shown in Eq. (1651) . We have also neglected all higher power corrections 

in q]_/P'j_ in the above derivations. From this result, we find the associated twist-three 

quark-gluon correlation function will take the form Tp{x,x — Xg) = Tf{x,^x), that is, its 

two arguments are not at the same. This is different from the soft pole contributions where 

Xi and X2 are equal in the twist-3 quark-gluon correlation matrix element. 

In our calculations, we further choose physical polarization for the radiated gluon, with 

the following polarization tensor: 

Y: e',{k')e,{kY* = -9'^' + ^p"^^ . (78) 

With this choice, we only have contributions from the three-gluon vertex diagrams 02, 62 
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FIG. 21: Illustration of the factorization of the hard pole contributions when k' is parallel to Pa- 
This part can he factorized into a hard pole contribution to the Sivers function in the polarized 
hadron A. The little circle in the diagrams represents all possible gluon attachments, as required by 
gauge-invariance. For example, the first one in the left panel represents all three diagrams 01-03 
of Fig. m 



and C2 and, apart from the color factor, each of them contributes a term j^. Their color 
factors are straightforward to decompose: 



1 1 



as : — — i-tfademT'^T''T''T'=)TT{T'T'T'^ 



1 1 



62 : i- — {-tfade)TT{T''T''T''rr)TT{T'r' 



1 1 
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2Nr 
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2Nr 



X C 



F2 



(79) 



We then directly find that the contributions of these diagrams can be factorized into the 
hard pole contributions to the Sivers function, multiplied by C/ (or Cfi and Cf2) and the 
partonic hard-scattering function. We display this factorization graphically in Fig. [2TJ The 
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total contribution is given by 

\hard 



af3 qa J3 



k'ocPA 



X 



271 



qq'^qq' (g^)2 

- (^-^ + Cf^ j dxi Tf{x, X-Xg 



1-e 



(80) 



where Xg = (1 — C,)x and H^^J'^g^, is given in Eq. ( 128|) . Indeed, the above result reproduces 
the contribution to Eq. fl38l) that arises from the hard pole part of the quark Sivers function. 

In summary, by combining all the soft-pole and hard-pole contributions when the radiated 
gluon is nearly parallel to the polarized nucleon, we can reproduce the contribution from the 
perturbatively generated Sivers function to the SSA in dijet-correlations in Eq. (j38ll . derived 
from the TMD factorization formula in the regime Aqcd q± P±- With these explicit 
calculations, we have demonstrated that at this perturbative order the contribution to the 
SSA when k' is parallel to Pa can be factorized into the quark Sivers function, multiplied 
by the initial/final state interaction color-factors (C/, Cpi, and Cf2), and the partonic 
scattering amplitude squared. 



3. Non-derivative terms: k' parallel to Pb 

The calculation for the contributions when k' is parallel to Pb follows the same procedure 
as above. As we have shown earlier, there is no contribution to the derivative terms when 
k' is parallel to Pb, and so we only need to consider the non-derivative contributions. We 
will also have both soft and hard poles contributions. In Fig. [151 we have shown the soft 
pole diagrams, and in Figs. [T7IIT8] the hard pole ones. 

Again we have to keep track of the fc^x-fiow in the scattering amplitudes. First, we need 
to parameterize the momenta of all the external particles including their linear dependence 
on kgj^. For example, when the polarized gluon attaches to the left of the cut line in the 
diagrams of Figs. [TSlfTSl the momentum /c^ of the radiated gluon and the quark momentum 
/c2 can be parameterized as 

kt = (1 - Pi + f^,^ii'fj^:p/ ^ + *!?x - . 
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Unp. Quark Distr. 



FIG. 22: Factorization of the contributions when k' is parallel to Pb- This contribution is factorized 
into the TMD quark distribution in the unpolarized nucleon B. The little circle in the diagrams 
represents all possible gluon attachments. The bars indicate the poles. The second and third dia- 
grams in the left panel represent the soft pole diagrams bi -65 and c\ -C5 of Fig. [73 respectively. The 
first one contains both soft poles (represented by the bar on the circle (diagrams 01-05 of Fig. [T^) 
and hard poles (indicated by the bar on the straight line (all diagrams in Fig. \1'T^). 

and ki = k2 — kg with kg determined by the pole. Note that these parameterizations are 
different from those in Eqs. fl73l) . fl7^ . This is because, when k' is parallel to Pb, its minus 
component is dominant and is parameterized by (1 —^')x'P^. Again, this minus component 
does not depend linearly on kg±. The plus component of k'j^ does have linear dependence, 
which can be calculated from the on-shell condition [k'^Y = 0. The linear /c^x-dependence 
for the plus component of can be calculated as well. 

Likewise, when the gluon attaches to the right of the cut, we will have 



and ^2 = ki + kg with kg determined from the position of the poles. Again, different poles will 
give different results for kg, following the above discussions for the soft pole contributions. 
By using the above parameterizations for the momenta of the external particles, we can 



2fci± -qi-qi- q± 
2x'{l - ^')Pb ■ Pa 




(82) 
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calculate the contributions from the soft pole diagrams in Fig. [151 the hard pole diagrams 
in Figs. [T7] and [181 We further find that the contribution from the hard pole diagrams in 
Fig. [TS]is power suppressed by q±/P±- So, for the hard pole contributions, we only need to 
consider the diagrams shown in Fig. [T71 which all have the same hard pole indicated by the 
bars on the propagators. As discussed in the last subsection, the hard pole contributions 
may in general lead to the polarized gluon having nonzero longitudinal momentum fraction, 
Xg 7^ 0. However, when the radiated gluon's momentum k' is nearly parallel to Pb, the hard 
poles lead to = at leading power, as was the case for the soft pole contributions. This 
is very different from what we found for the hard pole contributions in the last subsection 
in the case k' parallel to Pa- The momentum giving the hard pole, x'Pb — k' + kg with 
kg = XgPA, is the same as in the last subsection. However, taking the pole, we now find 

- 2k, ■ Pa ' ^ ^ 

where the last equation holds in the limit q±/P± — * 0. The reason for this is that, when k' 
becomes parallel to Pb, it will be dominated by its minus component and can be parameter- 
ized as k'^' = (1 - C)x'Pb + ^^a/2(1 - C)x'Pb ■ Pa + k'^. Substituting this into the above 
equation f[83|) . we find that the numerator is of order g^, whereas the denominator is of order 
P^. So, in the leading power approximation, ~ for the hard pole contributions when k' 
is nearly parallel to Pb- The contributions from the soft poles in Fig. [TB] and the hard poles 
in Fig. [T7] is added to obtain the final results. We emphasize again that our power-counting 
based analysis is only valid for the leading power terms in the q±/ Pi_ expansion. 

The final result for the contributions from Figs. [151 and [T71 plus their mirror diagrams, is 



dAer(5'_L)(ggO 
dyidy2dP]_d?qA 



_ _^ Sivers if, S -A ^ '^±q± 

k'ocPB 



Hq^,f'Z\q'{.^, t, u) ^2 '^Pi^a, Xa)] 



X ^C^l dx'^q'ix') 



1 + 



which is again proportional to the same hard factor H^^J'^gg,, and reproduces the last term 
in the TMD factorization formula Eq. (l38l) . It is important to note that a part of the 
hard pole contributions from Fig. [T7] cancels out the soft pole contributions from diagrams 
Oi-Os of Fig. [151 The remainder of the hard pole contribution is factorized into the quark 
distribution multiplied by Ci and the partonic scattering function. The soft pole diagrams 
61-65 and C1-C5 are separately factorized into the same quark distribution for the unpolarized 
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nucleoli multiplied by Cpi and Cf2, respectively, and the partonic scattering function. We 
illustrate this factorization in Fig. [22j 

After summing up all contributions from Eqs. (|7T]) . (175]) . (jHOD, and we obtain the 
total leading-power contribution to Aa{S±) from the {g)qq' qq'g partonic subprocess in 
the q±/P± expansion. It reproduces the factorized formula in Eq. (155]) . which is the leading 
term in the corresponding TMD factorization formula, Eq. ([2]). 

In summary, we have calculated the single transverse-spin asymmetry in dijet-correlation 
in the twist-three approach when the radiated gluon is parallel to either the polarized nucleon 
or the unpolarized nucleon. By using the power counting technique, we have shown that the 
contribution to the SSA can be factorized into the perturbatively generated quark Sivers 
function when the radiated gluon is parallel to the polarized nucleon, and into the unpolar- 
ized quark distribution when it is parallel to the unpolarized nucleon. We have demonstrated 
that the result reproduces the leading order terms in the corresponding TMD factorization 
formalism in the regime Aqcd q± <^ -Pl. This is a nontrivial result, especially for the 
single spin asymmetry, because the calculations for the quark Sivers function at this order 
have to take into account the perturbative expansion of the gauge link at next-to-leading 
order {0{g^)). 



IV. THE q2_ MOMENTS OF THE SSA 

The single spin dependent differential cross section term in Eq. ([2]) can be further sim- 
plified by taking a moment in q_i_. Such a moment was also considered in 22|. We can 
rewrite 

e'^^SIql = \S±\^-f^ sin (j)h sin 6 , (85) 

where (pb is the so-called bi-sector angle of the two jets: (pb = (0i + 02)/2 with 0i and 
02 the azimuthal angles of the two jets relative to the polarization vector S±. The angle 
6 = TT — (02"~0i) measures how far the two jets are away from the back-to-back configuration. 
All these azimuthal angles are defined in a frame in which the polarized proton is moving 
in the +z direction. 

Since |g_L| ~ |Pj_||sin5|, the g^-moment of the asymmetry is also related to the sin5- 
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moment. One interesting moment is the following: 



\P± 



■ sin 6 



ab 



-9 
Mp 



XamXa,Xa)x,Mx,)H!i;Zl',{Pl) , (86) 



where Tf{x,x) is the twist-3 matrix element of the quark-gluon correlation function defined 
in Eq. §B 



Eq. m 
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251. is ^Is*^ related to the fc_L-moment of the TMD quark Sivers function, see 
42|. An advantage of taking this moment of the asymmetry is that the trans- 



verse momentum integrals of the various factors in the factorization formula decouple from 
each other, without further assumptions for the fc^-dependence of the TMD distributions. 
In fact, the SSA in the q± moment of the dijet momentum imbalance is effectively a physical 
quantity with only a single large observed scale, P±. For such an observable, the conventional 
collinear factorization is more appropriate. In this case, all partons' transverse momentum 
dependence is integrated into the usual collinear parton distributions. The hard factors we 
obtained above for all the partonic channels agree with those given in 22], although a very 
different approach was adopted there. In 22(], in order to investigate the g_|_-moment of 
the single spin asymmetry, the appropriate gauge links were derived and expanded to first 
order in g. Since both approaches have included the initial/final state interactions in their 
formalisms, they should agree with each other at this order. 



V. CONCLUSION 



In this paper, we have studied the asymmetric production of two jets in hadronic colli- 
sions. Using the collinear factorization approach in perturbative QCD, we calculated both 
the spin-averaged and the spin-dependent differential cross section for the hadronic produc- 
tion of dijets with momentum imbalance q±, in a kinematic region where P± ^ q± ^ Aqcd- 
At the first nonvanishing order, the momentum imbalance is generated by radiating a gluon 
with transverse momentum equal to the imbalance. In the limit when the imbalance q± is 
much less than the averaged jet momentum Pi, we derived the leading-power contribution 
to the cross sections in the q±/P± expansion when the radiated gluon is nearly parallel to 
one of the incoming hadrons. We found that the perturbatively calculated leading contri- 
butions to both the spin-averaged and the spin-dependent cross sections can be factorized 
into a partonic hard part which is a function of Pi, and into perturbatively generated TMD 
parton distributions at 0{g^). Our results derived in the collinear factorization approach 
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in the limit q^/Pi_ — > reproduce the same factorized expressions in Eqs. (135|1 . (138!) that 
were derived as the leading order terms of the generalized TMD factorization formulas in 
Eqs. when the observed is solely due to the perturbatively generated TMD parton 

distribution of one of the incoming hadrons. 

The consistency between the results derived in the coUinear factorization approach and 
those derived by expanding the generalized TMD factorization formulas in Eqs. ([2]), ([5]) to 
the same order naturally leads us to ask if the generalized TMD factorization formulas in 
Eqs. ([2]), ([5]) are actually valid for the hadronic dijet production to all orders in perturbative 
QCD when Pi ^ gj. > Aqcd, like the TMD factorization formulas for Drell-Yan and SIDIS. 
// this were true, the consistency that we demonstrated in this paper in the overlap region 
where P± ^ q± would effectively provide a unified picture and treatment for the hadronic 
dijet production over the full kinematic range of jet momenta: using the TMD factorization 
formalism when the dijet momentum imbalance is small and the collinear factorization ap- 
proach when it is large. This is exactly what was achieved in Ref. 28| for the Drell-Yan and 
SIDIS processes. 

However, as we mentioned at the beginning of Sec. Ill, showing the consistency between 
two factorized formalisms at the leading order when q±/P± —>■ is not sufficient to prove 
the TMD factorization formula over its full kinematic region. When we started with the 
collinear factorization approach to calculate the momentum imbalance of the dijet cross 
section, we already took as a fact that the long-distance interaction between the hadrons/jets 
at their mass scales, O(Aqcd), are factorized for the given q± ^ P±. That is, our results 
derived in the limit q±/P± are valid when q± ^ Aqcd; but, our derivation does not 
address the issues concerning the factorization of the interactions between the hadrons/jets 
when g_L — > Aqcd or the mass scale of the hadrons/jets involved. The consistency that we 
demonstrated here is certainly a necessary condition that needs to be satisfied should the 
TMD factorization formulas be valid. 

For the TMD factorization to be true, we need to establish two key facts: the partonic 
hard part at 0{P±) must be insensitive to any physics at the scale q± ^ Pi, and the TMD 
parton distributions at the scale q± shoulc 
"quasi" -universal up to a sign change [10 
the partonic hard parts at 0{P±) by using the Brodsky-Hwang-Schmidt model with the 
proper color factor for the interaction vertices. We derived the partonic hard parts from 



be process independent and universal, or at least, 
In this paper, we also calculated in Sec. II 
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the lowest order scattering diagrams and showed that the dynamics of partonic scattering 
at scale 0{P±) is independent of the model of nucleon that we used in our derivation. Since 
the SSA can be generated by either initial- or final-state interactions, the color flow of the 
partonic scattering in the dijet momentum imbalance is very different from that of Drell-Yan 
and SIDIS. Since all poles from the initial- and final-state one gluon interactions contribute 
to the SSA, and the color factor for each individual scattering diagram is insensitive to the 
long-distance details of the nucleon because there is only one unique color structure for the 
active quark-gluon combination at this order, we were able to factor the process-dependent 
as well as the diagram-specific color factors along with the partonic hard parts, and leave 
the physics at scale 0{q±) in the TMD parton distributions defined in SIDIS. By doing 
that, we found that the resulting full partonic hard parts arc the same as those derived in 
the coUinear factorization approach. However, for a full TMD factorization, we will have to 
show that the poles from multiple gluon interactions can be exponentiated into the gauge 
link that defines the TMD parton distributions, while leaving the same partonic hard parts 
with the same color factors. Due to the mixture of initial- and final-state interactions and 
the complications of color fiows, this task is very non-trivial and is beyond the scope of this 
paper. 

To better understand the dijet momentum imbalance in hadronic collisions, a number 
of extensions can be performed based on our results. First, in this paper, we have only 
studied the contribution of one gluon radiation coUinear to one of the incident hadrons. 
This should be extended to other kinematic regions of the gluon. For example, when the 
gluon is radiated parallel to the final state jets or hadrons, it is important to see how the 
corresponding contribution can be accommodated by the jet definition, or be factorized into 
the fragmentation functions for the hadron, especially in the spin- dependent case. Another 
important kinematic regime is when the radiated gluon becomes soft. It is very important to 
see if the resulting contribution can be factorized into a soft factor, as predicted by the TMD 
factorization formalism. In addition, in this paper we have only discussed the SSA for dijet 
correlations, which has Sivers-type contributions. For the related di-hadron correlations, one 
expects also the Collins mechanism to be important, and it will be interesting to explore 
it and the factorization properties of the corresponding spin-dependent cross section. The 
method and procedure used in this paper should provide useful guidance in all these related 
studies. 
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We finally note that if the TMD factorization indeed fails, a careful comparison between 
our results and experimental data as the momentum imbalance q±/P± decreases could fur- 
nish an ideal test of the breaking of factorization and provide new opportunities to explore 
QCD dynamics. 
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APPENDIX A 



In this appendix, we list the hard factors for the dijet-correlation at hadron colliders, for 
the unpolarized and the single-transverse spin dependent cases. For unpolarized scattering, 
the hard factors are well-known in the literature. They may also be obtained from the results 
shown in Tables I-III, using Eq. ( 129|) . For completeness, we list them here: 
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For the single-transverse spin dependent case, using Eq. (l30|) . we obtain 
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The above results are for t-channel scattering. In dijet-production, both t and u channel 
processes will contribute. The hard factors for the -u-channel scattering can be obtained 
from the results above by exchanging t ^ u. 
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